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Theory of Riemann extensions of the spaces with constant afhne connection is proposed to study of 
the properties of nonhnear the first order systems of differential equations. 
As example quadratic system of differential equations 

where 

P{x, y) = ao + aix + a2y + anx^ + ai2xy + 022?;^, 

y) = 60 + bix + b2y + bux'^ + b^xy + 622?/^, 

and a, b are the parameters, is presented in homogeneous form and is considered as geodesic of three- 
dimensional space with constant affine connection depending on the parameters a, b. After the Riemann 
extension one get a six-dimensional space and its properties in relation to the parameters are investi- 
gated. 

The Lorenz system of equations 

dx , dy dz 

3- = cr(2/ a;), —^rx~y~xz, — ^ xy - bz 

ds as ds 

after presentation in homogeneous form is considered as geodesic equations of four dimensional space 
with constant affine connection. Based on the eight -dimensional Riemann extension of a given type 
space the properties of the Lorenz system are studied. 



1 Introduction 

The subject of consideration is the first order polynomial systems of differential equations 

= c'' + a)x^ + b'kx'x'' (2) 
as ■' ■' 

depending on the parameters a, 6, c. 

They play an important role in various branches of modern mathematics and its applications. 
In particular case of the system of two equations 

4^ = aO + alx + a2y + allx'^ + al2xy + a22y^, 
ds 

^ =bO + blx + b2y + bWx^ + bUxy + b22y^ 
ds 



2 FROM THE FIRST ORDER SYSTEM OF EQUATIONS TO THE SECOND ORDER SYSTEMS OF 0DE2 



there are many unsolved problems. 

The spatial first order system of differential equations 

dx , \ dy . dz \ /,\ 

— =P{x,y,z), —=Q{x,y,z), — = R{x,y,z) (3) 

with the functions P,Q,R polynomial on variables x,y,z are still more complicated object for the 
studying of their properties. 

As example the studying of the Lorenz system of equations 

dx dy dz 

-^^=<y-x), -=rx-y~xz, = xy - bz (4) 



and the Rossler system 



dx dy dz , 

— = — y — z, — = X + ay, — = ox — cz + xz (5) 
ds as ds 



which are the simplest examples of the spatial systems have chaotic behaviour at some values of pa- 
rameters represent the difficult task. 



2 From the first order system of equations to the second order 
systems of ODE 

The systems of the first order differential equations are not suitable object of consideration from the 
usually point of Ricmann geometry. 

The systems of the second order differential equations in form 

d'^x^ : , , dx^ dx^ 

are best suited to do that. 

They can be considered as geodesies of the affinely connected space in local coordinates x^ . The 
values n^j, = arc the coefficients of afhnc connections on . 

With the help of such coefficients can be constructed curvature tensor and others geometrical objects 
defined on variety . 

There are many possibilities to present a given system of the first order of equations in the form of 

(6). 

One of them is a following naive presentation. 
For the system 

f =P(x,i,), f =«(x,„) (7) 

after differentiation with respect the parameter (s) we get the second order system of differential equa- 
tions of the form (6) 

(P'x 1 dx dy dx 
'd^^P^ '^'ds^ UsUs 



d^y 1 dx dy dy 
d^^Q^^Us^^^^TsUs- 

Such type of the system contains the integral curves of the system (7) as part of its solutions and 
can be considered as the equations of geodesies of two dimensional space M'^{x, y) equipment by affine 
connections with coefficients 
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It is apparent that the properties of the system (7) have an influence on geometry of the variety 
M'^{x,y). 

Remark that a given system is equivalent the second order differential equation 
which has the solution in form 

dx ~ P' 

By analogy the spatial system of the first order differential equations can be written. 
As example the Lorenz system of equations is equivalent the second order ODE 



da;^ y \ax J \ X J ax x 
where 

b + a+1 ^ 1 b{a+l) , ^ + 1 b{r - 1) 

which can be obtained by the elimination of variable z from the system 

dy rx — y — xz dz xy — hz 

dx a{y — x) ' dx a{y — x) 

The second order ODE's 



^ + ai(x,y) j +3a2{x,y) +3a3{x,y)— + a4{x,y) = 0, 

with arbitrary coefficients ai{x, y) are form- invariant under the change of the coordinates 



and are equivalent to the system 

d'^x ^ ( dx\^ ^ 1^ ^ ( dx\ ( dy\ ^ ^ ( dy\^ ^ 



ds^ ' \ds I ' \ds J \ds J ' V , 

d'^y . . f dx\^ . ^ f dx\ f dy\ . . f dy\^ 

^+a4(..,y)^-j +2a3(x,,)(^-j(^-j+a2(:r,,)^-j =0, 

having the form of geodesies of two dimensional aSinely connected space ( with this aim the formulae 

d'^y xy — yx 



y = 


h{u, 


dx\ 


(dy 


Ts) 


\ds 


dx\ 


(dy 


d^j 


\ ds 



-) 

dx^ ' 



was used. 



3 From the afRnely connected space to the Riemann space 

Now we shall construct the Riemann space starting from a given aflinely connected space defined by 
the second order ODE's. 

With this aim we use the notion of the Riemann extension of nonriemannian space which was used 
earlier in the articles of author. 

Remind the basic properties of this construction. 
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With help of the coefficients of affine connection of a given n-dimensional space can be introduced 
2n-diniensional Riemann space D^" in local coordinates (a;% ^t) having the metric of form 



2"ds2 ^ _2nfe. [x'-)^kdx'dx^ + 2d^kdx'' (9) 

where are an additional coordinates. 

The important property of such type metric is that the geodesic equations of metric (9) decomposes 
into two parts 

x^ + U^^x'x^ = 0, (10) 

and 

-^+R[..x^x'^i = 0, (11) 

where 

ds ~ ds ' ds 

and are the curvature tensor of n-dimensional space with a given affine connection. 

The first part (10) of the full system is the system of equations for geodesic of basic space with local 
coordinates a;* and it do not contains the supplementary coordinates '^k- 

The second part (11) of the system has the form of linear N x N matrix system of second order 
ODE's for supplementary coordinates 

-^ + Ais)-+Bis)^ = 0. (12) 
Remark that the full system of geodesies has the first integral 

which is equivalent to the relation 



2*fc— = i/s + /x (14) 
as 

where /U, v are parameters. 

The geometry of extended space connects with geometry of basic space. For example the property 
of the space to be Ricci-flat i?^ = or symmetrical Riju-m = keeps also for the extended space. 

It is important to note that for extended space having the metric (9) all scalar curvature invariants 
are vanished. 

As consequence the properties of linear system of equation (11-12) depending from the the invariants 
oi N X N matrix-function 

„ „ IdA 1 
E = B - -— - -A^ 
2 ds 4 

under change of the coordinates can be of used for that. 

First applications the notion of extended spaces for the studying of nonlinear second order ODE's 
connected with nonlinear dynamical systems have been considered by author (V.Dryuma 2000-2008). 

4 Rigorous approach to geometry of planar systems 

The system of paths of two-dimensional space 5*2 in general form looks as 

x + Ul,ixf + 2Ul^xy + Ul^{yf=0, 
y + Uf^ixf + 2Ill^xy + Tll^{yf = 0, 

where the coefficients II^j = Uj^ . 
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The Riemann extension of the space ^2 is determined by the metric 

^ds'^ = -2zHl^dx^ - 2ml^dx^ - Azlll2dxdy - Atnf^dxdy- 

-2zll\2dy'^ - 2tn^2<^y^ + '2'dxdz + 2dydt. 

We shall apply such system of equations for the studying of the first order planar system of ODE's. 
With this aim we use the facts about the geometries of paths for which the equations of the paths 
admits independent linear first integrals. 
L.P.Eisenhart, 1925 

A necessary condition to geodesic admit the linear first integral 

, , dx'' 
0'i{x,y)—r- = const 
ds 

is 

where 



~ dxi 



and are the Christoffel symbols of the metric. 
We apply this conditions and their consequence 

where i?^^.; is the curvature tensor of the space to determination of the coefficients of equations nf^ 
using the vector ai in form 

a, = [Q{x,y),-P{x,y),Q,% 
By this means that the first order of equation 

dy_ ^ Q{x,y) 
dx P{x, y) 

or 

Q{x, y)dx - P{x, y)dy = 

is an integral of the paths equations. 

With the help of these conditions it is possible to state only three coefficients of afJine connections 



k 

For determination of others coefficients we use yet another the first order equation 

dy y{y - 1) 



dx x{x — 1) ' 

with a first integral 

yix) 

In this case a second vector bi is in form 



Cjx - 1) 
x-C ■ 



h = [yiy-l),x{x-l),0,0]. 

The equation plays an important role in theory of of planar the first order system of equations and 
was used in a famous article of Petrovsky-Landis (1956) ([?]). 
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So, from the conditions on the metrics to admits two hnear first integrals the coefficients n^^- of the 
paths equation are uniquely determined and have the form 

jji = {^Q{x,y))x{x-1) 

y^P{x, y) - yP{x, y) + Q{x, y)x'^ - Q{x, y)x 



(^P(x,y))x(a;-1) 



y'^P{x, y) - yP{x, y) + Q{x, y)x^ - Q{x, y)x 

and corresponding expressions for the coefficients Tl\i,\i\2,^22 

Remark that last two equations are reduced at the independent equations 

^ + M(s)^ + N{s)z{s) + F{s) = 

and 

^ + U{s)^^+V{s)t{s) + H{s) = 

with the help of the first integral of geodesies 

, ,dx , , dy .s 
'^'^Ts+'^'^d-s-^2-^ = ' 

of the metric (18). 

It is interested to note that such type of non homogeneous linear second order ODE's are connected 
with theory of first order systems of ODE's as the equations on the periods of corresponding Abel 
integrals. 

5 The second order ODE's cubic on the first derivative 

The first two equations of geodesic of the metric are equivalent to the one second order differential 
equation 



d^ {[^P 



3 



^ ^ i^^^y))^'-{lyPix,y))x){£y{x)) 
dx^ y^P{x, y) - yP{x, y) + Q{x, y)x'^ - Q{x, y)x 



1 



{{P:,-Qy)x^ + {Qy-P:,-2P)x + {Py)y' + {-2P-Py)y+2P) {£y^ 
y^P-yP+Qx^-Qx 

{-iQx)x^ + {Q.+2Q)x+{P,-Qy)y^ + {2Q-P, + Qy)y-2Q) £y 

y2p -yP + Qx^ - Qx 

-{lM^,y))y^ + {i,Q{x.y))y ^^^^ 



y^P{x, y) - yP{x, y) + Q{x, y)x'^ - Q{x, y)x 
The equation (14) is of the form 

d'^y f dy\'^ f dy\^ , dy 



dx^+^^^^'^Hd^J +3a3(x,y)-+a4(x,y)=0, 

and it has the invariants depending from the coefficients ai{x,y) (R.Liouville, 1880, T.Tresse, 1886) 
under transformations of variables {x, y). 

As it was shown by (E.Cartan, 1924) these invariants in general case are same with the invari- 
ants of two-dimensional surface V'^{x,y) in a four- dimensional projective space i?P^(^') (under the 
reparametrization and the change of coordinates ^*). 
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For example the invariants of equations with condition on coefficients 1^5 = same with the invariants 
of developing surfaces in a three-dimensional projective space 

The invariants of R.Liouville have been used successful in theory of ODE's and their applications in 
works author (V.Dryuma, 1984-). In particular the second order ODE with a Painleve property have 
the condition 2/5 = 0. 

In our case the equation has the particular integral 

dy ^ Q{x,y) 
dx P{x, y) 

and the function 

C{x - 1) 

as the first integral. 

From geometrical point of view the equation (14) corresponds two-dimensional surface in the i?P^- 
space. 

In this context it is interested to note the relation with the Petrovsky-Landis theory of limit cycles 
of the equation 

dy _ao + aix + a2y + anx'^ + ai^xy + a^iy^ 
dx bo + bix + 622/ + biix"^ + bi2xy + 622?/^ 

In the famous article oi P — L was developed approach to the studying of the problem of the limit 
cycles of the first order quadratic equation. 

Let us recall the basic facts of the Petrovsky-Landis theory. 

For studying of the closed curves of quadratic the first order equation is considered the equation 

^ ^ y{y - 1) 

dx x{x — 1) 

with solution 

> {x-C)- 

As it was showed for the closed curves the parameter C satisfies the algebraic equations 



5^a„(/ii)C" = 



where the coefficients an{lJ-i) are dependent from the parameters of the quadratic equation. 
Such type of equation arises from the condition 



{x - Cf [x{x - l)P(x, y) + y{y -l)Q{x,y)] ^^ ^ ^ 



x^{x — 1)^ 

which lead to determination of the closed integral curves of the first order system of equations 

dy _ ap + aix + a2y + aux^ + ai2xy + 022^^ 
dx bo + bix + 622/ + biix'^ + b^xy + 622?/^ 

As it was shown by Petrovsky-Landis in this case the parameter C is determined from calculations 
of the residues of integral after substitution of the expression 

_ C{x - 1) 

in it. 

The simplest of them are 
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{b2 + 2 bO + al2 + al + a2 + 2 aO + hl2 + hi ) (7^+ 

+ (-2 aO -hl2 - a2 - hi + 2 h22 -2hO)C -2h22 - h2 = 0, x = 0, 

{-h2-2b0- al2 - al - a2 - 2a0 - hl2 - hi) 
+ {a2 + hl + al2 + 2h0 + 2a0 - 2all)C + 2all + al =0, x=l, 

{hl2 - 2 h22) C +h2 +2h22 = 0, x = C. 



{-al2 + 2all)C -2all - al ^0, x^oo. 

Others equations on C from the combinations of these conditions are followed. 

According the results of (P — L) general quantity of the values C defined by a such type of equations 
is equal 14 and this number coincides with the quantity of closed solutions determined by the quadratic 
equation. 

On other side 11 curves from 14 can be transformed into the small neighborhood of the essential 
singular points of the equation 

dy_ ^ y{y- 1) 
dx x{x — 1) 

They are: (0, 1), (1, 0), (0, 0), (1, 1) and so on... 

As result only three closed curves do not be transformed into the neighborhood of the singular points 
and so the quantity of the limit cycles defined by the equation is equal three. 

It is interested to note that a some conditions on the parameter C are appeared in context of the 
second order ODE. 

In fact the result of joint consideration of the first order equations 

/ ^ Q{x,y) , ^ yjy- i) 

P{x,y)' x{x-l) 

the function 

C{x - 1) 

and the second order ODE (13) we get the conditions 

a{x, y)C' + p{x, y)C^ + ^{x, y)C^ + 6{x, y)C^ + e{x, y)C + ii{x, y) = 0, 

a{x,y) = {al2 + bl2)y'^+ 
+ {bl + b2 + {2a22 + 2all +2h22 + 2bll)x+ a2 + al)y+ 
+ {al2 + hl2) x^ + {al +b2 + bl + a2)x + 2a0 + 2 bO, 
(3{x, y) = 2 h22 y^ + ((-5 al2 -hl2 -2 h22) x + 2 b2 - hl2) 
+ (H h2-b al-3 a2-2 bll -6 b22-hl)x-a2) + 

+ (H hll+al2 -10 all-bl2-4 a22) x^ + 2 hO-2 b2-hl) 2/+ 
+(-2 al2+2 all)x^ + {-2 al2 -2 hi -2 al -3 bl2-2 a2) x'^-A b(M- 
+{-Sa0-2bl~Ah0~3b2~al-2a2)x-2a0, 
l_i{x,y) = -y^x^al2 + {-2x*^all + {al2 - al)x^)y + 2x^all +x^al. 



where 
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Prom these relation we get the relations between the parameter C and coefficients aij, bij, ai, bi 
As example at the values x = and a; = 1 we get conditions on the value C 

{a2 + al2 + bl + bl2 + b2 + 2 aO + al +2 bO) 

+ (-2 bO - 2a0 - bl - a2 -bl2 + 2 b22) C - b2 -2b22 = 0, 

{a2 + al2 + bl + bl2 + b2 + 2 aO + al + 2 bO) 
+ {2all -2bO - bl -2aO ~ al2 - a2)C -2 all -~ al = 0. 
The substitution {x = C\y = 1 — C) lead to the conditions on the value C 

{bl2 - 2 b22) C + 2b22 + b2 = 0. 

After substitution {y = 1 — x,x = 1,C = 1/Cl) we get 

{al + 2all) Cl'^ + {al2-2all + bl + a2 +2a0 + 2b0) Cl- 

-al - a2 - al2 - 2b0 - bl2 - bl - b2 - 2 aO = 0, 
and the substitution {y = 1 — x,x = 0,C = 1/Cl) load to the condition 

(-2 b22 - b2) + (-2 bO - bl2 - bl ~ a2 ~2aO + 2 b22) C1 + 

+al + a2 + al2 +2b0 + bl2 + bl + b2 + 2 aO ^ 0. 

In result of a such type consideration we have got the conditions of Petrovskii-Landis article on 
parameter C. 

6 Three dimensional homogeneous system 

As it was shown in article (V.Dryuma,2006) between the planar system 

^ = ao + aix + a2y + aux^ + ai2xy + a222/^, 
at 

dy 2 2 

— = bo + bix + b2y + bux + bi2xy + 6222/ , 
dt 

and a spatial homogeneous quadratic system of equations 

dx „, dy ^, , dz , 

— = P{x,y,z), — = Q{x,y,z), — = Q{x,y,z) 



of the form 



^ = 4oo-2^+4a2y2:+(3ai - b2)xz+4:a22y^ + {Sai2 - 2b22)xy+{2aii - bi2)x'^, 

^ = 460^^+461x2;+ (362 - ai)yz + 46iia;2 + (36i2 - 2aii)a;y+(2622 - al2)y^ 
dz 

— = -(ai + 62)2;^ - (2622 + 012)2/2; - (2aii + 6i2)a;2; 



exists some connections. 

For a such system the condition 



dP{x,y,z) ^ dQ{x,y,z) _^ dR{x,y,z) ^ ^ 
dx dy dz 
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is fulfilled and in variables 

t/.x x{t) y(t) 

it takes the form of the equation 

dri bo + biS, + 62?? + + ^12^?? + ^22?7^ 

equivalent the planar system . 

It is significant that the spatial system in the variables 

, , dX(t) , , dY{t) , , dZit) 

takes the form of geodesic equations of a three dimensional space 

d^X' „, dX^ dX^ 



,0 +'^\k—, = (15) 

dt^ dt dt ^ ' 

in local coordinates X' = iX{t), Y{t), Z{t). 

In doing so the coefficient T^-f. are constant and depend on the parameters a^, a^ , 6^, bij. 

The set of such coefficients can be considered as the coefficients of affinc connection on a three 
dimensional space in local coordinates X'. 

From geometrical point of view through the equations (5) on th(^ space the structure of affinely 
connected space with constant coefficients of connection F*^ is determined. 

7 Six-dimensional Riemann space 

For the n-dimensional space equipped with affine connection F^-^ the metrics of the Riemann extension 
V2n has the form 

2"ds2 ^ -2rikdX^dX''d^i - 2d^idX' (16) 

where Xi are an additional coordinates. 

Non zero coefficients of affine connections of the space H3 are 

rj3 = -4ao, r^3 = -2a2, PJa = ^(62 - 3ai), = -4a22, 



rl2 = ^(2&22 -3ai2), Fii = 6i2-2aii, 

1 



Vis = -Abo, rf3 = -26i, r23 = -(ai-362), Ffi = -46n, 

r'12 = 2^'^^ll ~ 3^'l2)) r22 = ^12 — 2622) 

= 01 + 62, = ^(ai2 + 2622), r?3 = ^{2an + 612). 
According with the (16) the Riemann metric of six dimensional extended space Vq has the form 

^ds^ = {-2Ubl2 + AUall +8bll V) dx^ + {8 a22 U - 2V al2 + AVb22) dy^+ 

+ {8bOV-2 Wb2 + 8aOU -2 Wal ) dz^ + (-4 Ub22 - 4 Vail + 6 Vbl2 + 6 Ual2) dx dy+ 
+ (6 Ual -2Ub2-A Wall - 2 Wbl2 + 8blV) dx dz+ 
+ (-2 Wal2 + 8a2U + 6 Vb2 -2Val -4 Wb22) dy dz + 2 dx dU + 2 dy dV + 2 dz dW. 
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Remark that we use the denotes (x, y, z) for the coordinates {X, Y, Z). 

Geodesic of the metric consist from two parts. NonUnear system of coupled equations on coordinates 

{x,y,z) 

^x{s) + (-2 all + bl2) (^^(s)) + 2 (-3/2 al2 + h22) (^^a;(s)^ 
+2 (-3/2 +1/2 6.) (^,(,)) j-^z[s)-,a22 {^yis))^ -^a2 (^±yis)) ^zis^AaO (^.(.^ ' 
^,(.) - 4 6ii (^.(,)) ' + 2 (-3/2 + all) ^^yis) - 4 M (^.(s)) 

+ (-2 h22 + al2) (^2/(s)) + 2 (-3/2 6. + 1/2 al) (^2/(s)) ^^(s) - 4 W (^^w) = 0- 

^^(5)+2 (aii + 1/2 hl2) {j;x{s)^ i (V2 + 6^^) ^^(^)+(«^ + «'^) (^-^(«) 

And the Hnear system of equations on coordinates {U, V, W) 

^U{s) + Al ^U{s) + Bl ^V{s) + CI ^W{s) + El U{s) + Fl V{s) + HI W{s) = 

Clo do (Jib do 

St^V(s) + A2 ^U(s) + B2 ^V(s) + C2 ^W(s) + E2 U(s) + F2 V(s) + H2 W(s) = 
cts^ as as as 

-^W(s) + A3 -^U(s) + B3 ^V(s) + C3 -^W(s) + E3 Uis) + F3 V(s) + H3 Wis) = 0, 

ds^ as as as 

where the cocfBcicnts (Ai, Bi, Ci) are depended from the parameters (a, 6) and derivatives (i, y, z) 
with respect to the parameter s. 

Taking in consideration x, y, 2;-equations the U, V, W- equations can be one time integrated and 
take the form 

-^U{s) = (2 all + bl2) z(s)W(s)+ 
as 

+ ((-3 bl2 + 2 all) y{s) - 8 bll a;(s) -4bl z{s)) y(s)+ 
+ ((2 bl2 -A all) x(s) + (-3 al2 + 2b22) y{s) + {b2-3al) z{s)) U{s), 

-^V(s) = (al2 + 2 b22) z(s)W(s)+ 
as 

+ ((-3 al2 + 2 b22) x{s) - 8 a22y{s) - 4 a2z{s)) U{s)+ 
+ ((-3 bl2 + 2 all) x{s) + (2 al2 - 4 b22) y{s) + {-3 b2 + al)z{s)) V{s), 

-^W(s) = ((b2-3al) x{s) - 8 aOz{s) - 4:a2 y{s)) U(s)+ 
as 

+ {{-3b2 + al)y{s) -8b0z{s) -4:blxis))V{s) + 

+ ((2 all + bl2) xis) + {al2 + 2 b22) y{s) + 2 (ai + b2) z{s)) W{s). 

In result we have got a six-dimensional Riemann space associated with a quadratic the first order 
system of equations. 

An investigation the properties of the metric at the change of parameters a, b may be useful for the 
theory of a such type of the systems. 

In particular a study of the Killing properties of the metric allow us to get information on particular 
integrals of geodesic equations. 
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8 On the surfaces defined by spatial system of equations 

The equation of surfaces z = z{x,y) defined by the spatial system of equations (4) has the form of the 
first order p.d.e 

(4 aO {zf + {4:a2y+ {3al - b2) x) z + A a22 y^^ + z^ ((3 al2 - 2 b22) xy + (2 all - hl2) x^) + 

+Zy (460^^ + {2,h2 - al)yz + A hll x^ +Ahl xz) + Zy ((2 b22 - al2) + (-2 all + 3 hl2) xy) + 

+ {h2 + al) {zf + {{2b22 + al2)y + {2all + hl2)x)z. (17) 

An examples of solutions of this equation were obtained by the method of (u, w)-transformation 

developed earlier by author. 

To integrate the partial nonlinear first order differential equation 

F{x,y,z{x,y),Zx,Zy) = Q (18) 

we use a following change of the functions and variables 

z{x,y) ^ u{x,t), y^v{x,t), z^^u^-^Ut, Vy^^. (19) 

In result instead of the equation (18) one get the relation between the new variables u{x, t) and v{x, t) 
and their partial derivatives 

^{u,v,u^,Ut,Vx,Vt) =0. (20) 

In many cases the solution of last equation is a more simple problem than solution of the equation 
(18). 

In result of application of the (u, w)-transformation the equation for z = z{x,y) takes the form of 
(20). 

Using the substitution 

d d 
u{x, t) = t-g^uj{x, t) - uj{x, t), v{x, t) = ^w(a;, t), 



where u{x, t) = xA{t), we find from a given relation the equation for the function A{t) 
{-t^h2 -th22 + a22 A{t) -t^hO + aO A{t)t^ + a2tA{t)) (^^^W^ + 
+ (-2 aO {A{t)f t - tbl2 - a2 {A{t)f -t^bl + tal A{t)^ ^^(*)+ 



+ (2 t^bO A{t) + m A{t) + al2 A{t)) ^A{t)- 



-tbll-tbO iA{t)y-al iA{t)y+tblA{t) + allA{t) + aO {A{t))=0. 

Solutions of this equation depend from the parameters a, b and can play a key value in theory of 
the quadratic systems. 

Another type of substitution 

d d 
v{x, t) = t-g^<^{x, t) - uj{x, t), u{x, t) = ^w(a;, <), 

where u){x,t) = xB{t), lead to the equation on the function B{t) 

{a22 B{t)t^ + a2 B{t)t + aO B{t) + b22 + tb2 + bO) { j^B{t) 
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+ (tbl2 - a2 {B{t)f + bl +al B{t) - 2 b22 tB{t) - 2 a22 {B{t)f t + tal2 B{t) - b2 B{t)j j^B{t)- 

-al2 {B{t)f - bl2 B{t) + all B{t) + a22 {B{t)f + b22 {B{t)f + bll =0 (21) 

which is a more simple than previous and also can be useful in theory of quadratic systems. 
Let us consider some examples. 
In the variables 

the equation (21) determines algebraic curve 

F{S,T,t,a,b) = 

having genus g' = 1 or g' = in depending on the parameters a, b. 
As example at the conditions 

60 = 0, aO = 0, 62 = 0, al = 0, a2 = 0, a22 = 0, 611 = 0, al2 = 0, all = 612, 61 = 

the equation (21) takes the form 

b22 (^A{t)) + (-2 b22 tA{t) + tbl2) ^-A{t) + b22 {A{t)f = 

and defines algebraic curve of the genus ,9 = 0. 

Integral of this equation is defined the relation 



ln(i) + ln(-^-4 b22 B{t)bl2 + bl2'^ - bl2) 



bl2 



^J-A b22 B{t)bl2 + bl2^ - bl2 

I J) I ' hi 2 

- ln(^-4 b22 B{t)hl2 + bl2^ + bl2) - 1/2 — — - - ln(i3(i))- 

^-4b22B{t)bl2 + bl2^ + bl2 
-.CI = 0. 

After the inverse (u, v)- transformation with the help of the function B{t) can be find the solution 
of the equation (17) at indicated above the values of parameters. 

9 Eight-dimensional Riemann space for the Lorenz system of 
equations 

To investigation of the properties of classical Lorenz equations 

dx , . dy dz ^ 

— =a{y-x), — =rx-y-xz, — = -bz + xy (22) 
ds ds ds 

we use its presentation in the form 

^C(s) = 1/5 {b-4a+l)^p + ar]p, 

-^V{s) = + {b + a-4)r]p, 

ds 

-^9{s)=^r]+ 1/5 ia-Ab+l)p0, 
ds 
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^p(,s) = l/5 + l + 
The relation between both systems is defined by the conditions 

xis) = -, vis) = -, z{s) = -. 
P P P 

Four dimensional system can be presented in the form 

which allow us to consider it as geodesic equations of the space with constant affine connection. 
Nonzero components of connection are 

1 _ Aa-h-2 pi _ cr p2 _ 1 

il4- ^ , ^24 --2' ^13-2' 

_9 r _o 4 — (T — 6 _o Ah — a — 1 

il4- 2' ^34- , ^34- ' 

3 _ 1 4 _ a+b+1 

il2--2' -^44- ^ 

The metric of associated space is 

^ds^ 2/5 (fe + o- + 1) du^V+ 

+ (2/5 dzdu + 2 dx dy - 8/5 dz dub + 2/5 dz du a) U+ 
+ (2 rdx du + 2/5 dz du a — 2 dx dz — 8/5 rfz (iu + 2/5 dz du b) Q+ 
+ (-8/5 dx dua + 2a dy du + 2/5 dxdub + 2/5 rfa; dw) P+ 

+2rf?/dO + 2(izrf{/ + 2rfMrfF + 2da;dP. (23) 
After integration geodesic of additional coordinates take the form 

^P(i) - {z{t) - ru{t)) Q{t) - 1 (-6 - 1 + 4a ) u{t)P{t) + y{t)U{t) = 0, 

^Q(i) + x{t)U{t) +au{t)P{t) = 0, 

iu{t) + i ((tr + 6 - 4) M(t) - xit)) Q{t) + \ (-4& + u{t)U{t) = 0, 
at b 5 

Ay(i) + 1/5 (6 + . + 1) uitmt) - --mmi + ^wQWfW ^ 0. 

Ctf U(f) '(i(f) 

Properties of this system from the parameters are dependent and can be investigated with the help 
of the Wilczynski invariants. 
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10 Laplace operator 

In theory of Riemann spaces the equation 

can be used to the study of the properties of the space. 

For the eight-dimensional space with metric (23) corresponded the Lorenz system we get the equation 
on the function 'tjj{x, y, z, u, P, Q, U, V) 

( S:r-_th \ P - 2a pS:r-_^ - 2 ( S:r-_tb\ rQ + 8/5 ( -rjL.^^ Q-2/5 ( 



-r^^ib ]Pb-2/5V -^ib 6 - 2/5 y -r^^ cr - 2/5 V-r^^b - 2/5 - 



-^/^ [dP9V^ ) - ""[W^V'- ^ lay^ ^ j - V5 - 2/5 J Q6- 

This equation has varies type of particular solutions. 
As example 

i;{P,Q,U,V) = H{P,Q,U) + VP 
where the function ipiPj U, V) satisfies the equation 

2 Qq^H{P, Q, U) + 9/5 Pa - 2 U^^H{P, Q, U) - 1/5 Pb-2rQ-6/5P = 0. 
Its solution is in form 

H{P, Q, U) = (^-3/10 - 1/20 b+^a^ arctan(|)p2 + PrU + _F1 (Q, U) + .F2{P, + t/^) 

where -Fi are arbitrary functions. 
At the condition 

6 = 9cr-6 

this solution takes a more simple form. 
In the case 

^{P,Q,U,V) = H{P,Q) + V 

we get the solution 

HiP. Q) = .F2iP) + .« ,0, + (lA°'-ys^+V'0<^)^Q 

for which the relation 

6 = 4-0- 

is a special. 
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11 Eikonal equation 

Solution of the eikonal equation 



dF dF 



also give us useful information about the properties of Riemann space. 
In the case of the space with the metric (23) we get the equation 



MIp') ij'^''' ifp') (^*) Hp*) H') Hp*) H 

- Hp*) H*)'''-'<-k*) w*-^^' H*) H*)^-'^' H*) H*)^' 
-/^ H*) Hv*) H*) H*) H*) H*) "'-'I' H 

-2/5 y ( A^) \ - 2/5 V ( A*) % - 2/5 v ( A 



one 



In particular case 

<i){x, y, z, u, P, Q, U, V) = A{Q, U, V) + PU 
get the equations to determination of the function A{Q, U, V) 



^A{Q, U, V)+2 QU-^A{Q, U, V)+2 QUP+6/5 U (^^(Q, U, V)^ Pa+6/5 U (^^(0, U, V)^ Pb- 

-4/5 U (^^(Q, U,V)^P-2U (^^(0, U, 1/)) rO - 2 aP (^^(Q, U, V)^ ^MQ, U, V) + 
+8/5 (^-^A{Q,U,V)^ (^A.AiQ,U,V)^Q + 8/5 (J^A{Q,U,V)^ QP- 
-2/5 (^-l^A{Q,U,V?j I^^A{Q,U,V)j Qa-2/5 (^-l^A{Q,U,V)^ QaP- 
-2/5 (^-^A{g,U,V)^ (^±.A{Q,U,V)^ g6 - 2/5 (^-^A{Q,U,V)^ QbP+ 
-8/5 (^±AiQ,U,V)^ (^A.A{Q,U,V)^ Ub-2/5 (^-^A{Q,U,V)^ (^l-AiQ,U,V)^ Ua~ 
-2/5 (^±A{Q,U,V)^ (^A.A{Q,U,V)^U-2/5v(^^A{Q,U,V)^ i,-2/5V (^-^AiQ,U,V)^ a- 

-2/5 y(AA(g, [7,^)^ = 

The solution of this equation is 

A{Q,U,V) = 

= .F00{U).C3V+l/2 2lnrrn +3/5 QU+3/5 ^-2/5 ^+2/5 ^-1/10 Q^-1/10 9^+_F5{U) 

-C^-rUU\U)a a G a o 
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where 



_FOO{U) = 5 



U 



-4_C3 + cr_C3 + _C36' 



.F5{U) = 1/2 




and 



6 = 2/3 + 2/3ct, r = 1/3+1/30-. 
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